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Abstract. In this article we argue that thermal reservoirs (baths) are 
potentially useful resources in processes involving atoms interacting with 
quantized electromagnetic fields and their applications to quantum technologies. 
One may try to suppress the bath effects by means of dynamical control, but 
such control does not always yield the desired results. We wish instead to take 
advantage of bath effects, that do not obliterate “quantumness” in the system- 
bath compound. To this end, three possible approaches have been pursued by 
us: (i) Control of a quantum system faster than the correlation time of the bath 
to which it couples: Such control allows us to reveal quasi-reversible/coherent 
dynamical phenomena of quantum open systems, manifest by the quantum Zeno 
or anti-Zeno effects (QZE or AZE, respectively). Dynamical control methods 
based on the QZE are aimed not only at protecting the quantumness of the 
system, but also diagnosing the bath spectra or transferring quantum information 
via noisy media. By contrast, AZE-based control is useful for fast cooling of 
thermalized quantum systems, (ii) Engineering the coupling of quantum systems 
to selected bath modes: This approach, based on field -atom coupling control in 
cavities, waveguides and photonic band structures, allows to drastically enhance 
the strength and range of atom-atom coupling through the mediation of the 
selected bath modes. More dramatically, it allows us to achieve bath-induced 
entanglement that may appear paradoxical if one takes the conventional view that 
coupling to baths destroys quantumness. (Hi) Engineering baths with appropriate 
non-flat spectra: This approach is a prerequisite for the construction of the 
simplest and most efficient quantum heat machines (engines and refrigerators). 
We may thus conclude that often thermal baths are “more friends than foes” in 
quantum technologies. 
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In this article we wish to ^^make out a case’’ for thermal reservoirs (baths) as 
potentially useful resources in quantum optics [1], namely, in processes involving 
matter interacting with quantized electromagnetic fields, and their applications to 
quantum technologies: quantum information processing [2, 3, 4], quantum sensing and 
metrology [5, 6, 7, 8, 9], as well as quantum thermodynamics [10, 11, 12]. In general, 
there is little we can do to avoid the ubiquitous presence of environments described 
as thermal baths in contact with quantum systems: with very few exceptions, all 
quantum systems are open [13, 14]. One may try to reduce the bath effects on the 
quantum system of interest by means of dynamical control, originally developed to 
suppress bath-induced decoherence or dissipation [15, 16, 17, 18, 19, 20, 21]. Yet such 
control does not always yield the desired results, hence we wish to advocate a different 
strategy that may be colloquially summarized as follows: Y/ you can’t fight the hath - 
join it”^ namely, take advantage of its effects, particularly those that do not obliterate 
“quantumness” in the system-bath compound. To this end, three possible approaches 
may be pursued, to be discussed in the subsequent sections p. 

• Control a quantum system faster than the correlation (memory) time 
of the bath to which it couples: Such control allows us to reveal quasi- 
reversible/ coherent dynamical phenomena of quantum open systems, manifest 
by the quantum Zeno or anti-Zeno effects (QZE or AZE, respectively) [28, 29, 
30, 31, 32, 33]. Dynamical control methods based on the QZE are aimed at 
protecting the quantumness of the system [30, 34, 35, 36, 37, 38, 39, 33, 40, 41, 
42, 43, 44, 45, 46, 47], but also diagnosing the bath spectra and transferring 
quantum information via noisy media (Sec. 2). By contrast, AZE-based control 
is useful for fast cooling of thermalized quantum systems [48, 49, 50] (Sec. 5). 

• Engineer the coupling of quantum systems to selected bath modes: 

This approach, based on field -atom coupling control in cavities [51, 52, 53] and 
photonic band structures [54, 55, 56, 57, 58, 59], allows to drastically modify 
bath-mediated exchange of virtual quanta between quantum systems and thereby 
enormously enhance their coupling [60, 61, 62]. Not less dramatically, such 
engineering allows us to achieve bath-induced entanglement [63, 64, 65, 66, 67, 
68, 69, 70] that may appear paradoxical if one takes the conventional view that 
coupling to baths destroys quantumness [13, 14] (Sec. 3-4). 

• Select or engineer baths with appropriate non-flat spectra: This approach 
is a prerequisite for the construction of the simplest and most efficient quantum 
heat machines (engines and refrigerators) [71, 72, 73] and for investigating their 
ability to attain the absolute zero [73] (Sec. 6). 

Our conclusions and outlook to forthcoming research along the discussed lines are 
presented in Sec. 7. 

2. Control within the bath memory-time: Zeno anti-Zeno dynamics 

Our theory of quantum systems whose weak interaction with thermal baths is 
dynamically controlled [30, 34, 40, 42, 44, 74, 75, 76] treats all kinds of such control, 

t Cooperative (Dicke) effects mediated by the bath are outside the scope of this article - cf. the 
following articles:[22, 23, 24, 25, 26, 27] 
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be it coherent or projective (non-unitary), continuous or pulsed, as generalized forms 
of two generic effects or control paradigms. One is 

Minimized hath effect = Quantum Zeno effect (QZE), 

which minimizes (under constraints on the control energy) the integral product 
(overlap) of two functions: G(co’), the coupling spectrum of the bath (obtained by 
Fourier-transforming its autocorrelation function) and a spectral “filter” function Ft(uj) 
determined by the control field-intensity spectrum and its time duration t. It is the 
“filter” function that provides the control handle on our ability to optimally execute a 
desired task in the presence of a given bath. In the presence of several baths (a common 
situation), both G{uj) and Ft{uj) functionals are represented by matrices [74, 75, 76]. 

QZE-based control is required in operational tasks related to quantum information 
its storage and transmission [77, 78, 79, 80, 81], where bath effects are detrimental 
and must be suppressed. Regardless of the chosen form of control, the controlled- 
system dynamics must then be Zeno-like^ namely, result in suppressed system-bath 
interaction. 

The alternative paradigm is 

Maximized hath effects Anti-Zeno effect (AZE), 

which amounts to maximized overlap of G{uj) and Ft{uj) (under control-energy 
constraints, as for QZE). AZE-based control is instrumental for non-unitary operations 
that entail changes of the system’s entropy. Such operations benefit from efficient 
interaction with a bath for their execution. Examples are measurements used to cool 
(purify) a quantum system [50], equilibrate (thermalize) it with a bath [49, 75], or 
harvest energy from the bath. If the underlying dynamics is anti-Zeno-like [29, 48], 
system-bath interaction will be enhanced and thereby facilitate these tasks. 

Certain tasks may involve state transfer or entanglement via the bath, which 
require maximized bipartite coupling, but minimized single-partite coupling with the 
bath [43, 74, 82, 83]. Eor such tasks, a more subtle interplay of Zeno and anti-Zeno 
dynamics may be optimal and depend on the quantum statistics of the bath [84]. 

We have therefore developed a general approach that allows to optimize the 
interaction of a quantum system with the environment so as to execute a given 
operation, be it non-unitary or unitary, such as state- transfer or storage with 
maximized fidelity, purification/entropy-minimization, entanglement distribution, 
or energy transfer [75, 76]. This approach consists in designing the temporal 
dependence of the Hamiltonian that governs the system by variational minimization or 
maximization (as the case may be) of a state-dependent functional chosen to quantify 
the success probability of the operation. To this end, the temporal control must be 
faster than the bath correlation time [38, 75, 76]. This approach not only provides 
protection from adverse effects of the bath, namely, quantum-state decoherence, but 
actually benefits from the system-bath interactions for the realization of a given non- 
unitary task. More formally, it maximizes the fidelity of any given quantum operation 
on a multidimensional Hilbert space for the baths or noise sources at hand. Its main 
merit is that it is not restricted to pulsed forms of control, and therefore can drastically 
reduce the energy required to execute a task by resorting to a smoothly varying field, 
thereby reducing the errors incurred by control [75]. 
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We shall discuss applications of dynamical control of the system-bath coupling that 
go beyond its conventional use as a means of fighting decoherence [15, 16, 17, 18, 
19, 20, 21, 34, 38, 39, 40, 85]. The first application of such control is as a tool of 
bath-spectrum diagnostics. Such diagnostics has the goal of revealing the dynamics 
of decoherence processes and their underlying bipartite and multipartite interactions 
(collisions). 

The basis for this diagnostics is the Kofman-Kurizki (KK) universal formula 
[29, 30, 34] 


R{t) = j Ft{uj)G{uj)duj change Ft (cj) ^ infer G(c(;) . (1) 

Filter Coupling—spectrum 

The diagnostic method consists in changing the filter function Ftiuj)^ e.g.^ by varying 
the control-field Rabi frequency, recording the resulting decoherence rate R{t) and 
deducing G{uj) from Eq. (1). To this end, the system, e.g. a qubit, is initially taken 
to be in a superposition of its excited (|t)) and ground (|^)) energy states. This initial 
superposition state 

|'0(O)) = cos 0 It) + sin 0 |t), (2) 

is subject to bath-induced decoherence (pure dephasing). It then has, at time t, a 
mean coherence that decays in a fashion dependent on R{t) 

(cra;(0)), (3) 

which is inferred from the probabilities of measuring the system in the symmetric or 
antisymmetric superpositions of energy states (Fig. 1) 


p±{t) 


1 ± 

2 


( 4 ) 


We have demonstrated (in collaboration with Davidson’s group) [86] the ability to 
infer the bath-coupling spectrum G{uj) via formula (1) by measurements performed on 
a large ensemble of cold atoms in an optical trap. A field with narrow spectral band was 
used to realize a filter function Ft{uj) that scanned the overlap integral in Eq. (1) upon 
varying the field strength (Rabi frequency). By measuring the decoherence rate R(t) 
as a function of the filter value we could infer the bath-coupling spectrum in the weak- 
coupling limit. This demonstration has experimentally established that the Kofman- 
Kurizki (KK) universal formula (1) allows the design of dynamical control (continuous- 
wave or pulse sequence) that is optimally adapted to the measured coupling spectrum 
of the bath. 


2.2. Maximized information on the hath by dynamical control 

We have recently been studying the maximum information obtainable on unknown 
spectral parameters of a bath (environment) by controlled spin qubits that serve as 
its probes [87]. This information is important for maximizing the sensitivity of spin 
probes at nanoscales, serving as magnetometers, thermometers, sensors for imaging 
or monitoring chemical and biological processes [88, 89, 90, 91]. 

By using tools of quantum estimation theory, we can find the precision of 
estimating key parameters of environmental noises (baths) that the spin (qubit) 
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Figure 1. Schematic view of probing bath parameters by a qubit that 
undergoes bath-induced decoherence while being subject to dynamical control. 
Information on the bath that yields estimation of its parameters is extracted 
from a measurement of the qubit in the |zb} = ^ l-^)) basis of Eq. (2). 


can probe. These include the probe-bath coupling strength the correlation time 
of generic bath spectra Tc, as well as their dephasing time T 2 . By optimizing the 
dynamical control on the probe under realistic constraints one may achieve the best 
accuracy of estimating these parameters by the least number of measurements possible. 

To this end, we minimize the relative error of estimating a bath parameter xb by 
means of the dependence of the decoherence rate of a qubit-probe R{t) = R{xB,t) on 
G{uj) = G{xb^oj)^ as described by Eqs. (l)-(4) (Fig. 1). This error obeys the bound 


, SXB . 1 

S[XB^t) = - > - 

Xb XB\/NmJ^Q(xB,t) 


( 5 ) 


Here we have introduced the number of measurements, and the quantum Fisher 
Information (QFI) [92, 93, 94, 95] for the qubit probe that is subject to dephasing as 
well as dynamical control 


^—2R(xB,t) t 

TQ{xB,t) = 


/ dR{xB,t)t 

V dxB 


2 


(6) 


where 0 is as in Eq. (2). In general, we can minimize the relative error per 
measurement by maximizing QFI: 


d^Q{xB:topt) = max?FQ{xB,t) 


( 7 ) 


which amounts to preparing the optimal initial state (2) with 0 = J, measuring the 
qubit at the optimal time and efficiently controlling the quantum probe. 

To demonstrate the potential of this approach, we may estimate, for example, the 
correlation time Tc, a key parameter of Ornstein-Uhlenbeck processes characterized by 
Lorenztian bath spectra. 


G{xb,(^) 


g^rc 


7r(l + ’ 


( 8 ) 


assuming that the system-bath coupling strength g is known. 

Dynamical control of the qubit probe can drastically improve the estimation of 
Tc'. as shown in Fig. 2, sequences of equidistant 7r-pulses (phase-flips), known as Carr- 
Purcell-Meiboom-Gill (CPMG) sequences [96, 97, 98, 99], give rise to a minimal error 
estimation that obeys ScPMci^c^topt) — while under free evolution the error 

V m 

grows with gxc^ Sfreei'T'c^topt) ^ 1 therefore much larger for gxc ^ 1, i.e. 

for distinctly non-Markovian bath spectra. 





















Thermal Baths as Quantum Resources: More Friends than Foes? 


6 



time t time t 


Figure 2. Estimation accuracy of the spectral width Tc^ of an Lorentzian bath 
with a quantum probe freely evolving (left) and dynamically controlled with a 
CPMG sequence oi N = ^ 7r-pulses (right). The minimal error per measurement 
^('rc, topt) is much smaller for a qubit probe under optimal dynamical control than 
under free evolution (middle). 


2.3. Bath-mediated transfer of quantum information 

The ability to transfer an unknown quantum state between nodes where the quantum 
information (QI) can be reliably stored and/ or processed is at the heart of QI 
processing and communication schemes. Since practically any medium connecting 
distant nodes corrupts the QI [100, 101, 102], one commonly resorts to probabilistic 
quantum repeaters [103], effected by conditional measurements [104]: only the desired 
outcomes are kept while the undesired outcomes are discarded. Such protocols [103] 
are severely constrained by high qubit-overhead and long average duration of successful 
QI transfer. It is clearly desirable to resort to deterministic protocols whenever 
possible. Here we advocate the possibility of such protocols, whose high success rate 
relies on dynamical control that is optimally adapted to the medium [79, 80, 81]. 

The idea is to write the full Hamiltonian as 


H = HsFHbFHsb. (9) 

Here the system S consists of the two spins that constitute the nodes between which 
the QI is transferred and a mode (channel) of the medium that couples these spins, 
labeled by /c = z, all other modes being treated as a thermal bath B to which S is 
coupled. 

The transfer fidelity over time T is again governed by the KK universal formula 
[34, 79] 

/ c» 

dujFT{oj)G{Lo). ( 10 ) 

-OO 

We need not know the detailed spectral distribution of the S — B coupling G{uj), 
only its width and crude mode spacing, which can be estimated by the methods of 
Secs. 2.1-2.2. Such estimation should suffice for designing the optimal tradeoff of the 
fidelity F versus time transfer T by appropriate temporal modulation of the coupling. 
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Figure 3. (Left) Schematic view of quantum information transfer through a 
fluctuating quantum spin-channel with random coupling J. The transfer over 
time T is optimized by an appropriate modulation Q:(t), |Q:(t)| < 1, of the 
boundary qubit-couplings to their nearest-neighbors. (Middle) Schematic view of 
the coupling spectrum between the “system” (comprised of the two boundary 

qubits 0 and N 1 and the spin-chain mode that couples them) and all other 
modes of the spin chain viewed as a bath (red) and the Alter Ft{oo) corresponding 
to the optimal chosen modulation shape a{t) oc (^) (green) with chopped off 
tails, as compared to the (free evolution) Alter (black), which has extended tails 
and therefore much larger overlap with G(uj). (Right) The chopping-off of the 
tails in the (green) Alter may reduce by several orders of magnitude the transfer 
infldelity (error) as well as the time required for transfer. 


Strikingly, one may analytically prove [79], upon parameterizing the modulation 

a{t) = aosinP p = 0,1,2 \ao\<l (11) 

that the best tradeoff is usually achievable for p = 2, because it yields a filter 
function without spectral tails (these vanish abruptly at a controllable frequency). 
This chopping-off of the tails may reduce by several orders of magnitude the transfer 
infidelity (error) as well as the time required for transfer (Fig. 3)! 

This control method, which is universally applicable to media consisting of 
interacting fermions (spin-| particles) and bosons alike, may also be used to maximize 
the storage time of QI inside a bath memory embodied by an inhomogeneously 
broadened and thermally fluctuating spin ensemble [81]. Thus, the coupling between 
quantum systems via the bath is required for effecting QI transfer or storage, and 
system-bath interaction control may serve as a tool for optimizing these processes, 
on the basis of minimal knowledge concerning the bath. This method is another 
application of our universal procedure for fidelity optimization of the task at hand 
and the ability to prioritize the use of resources for implementing it in any given bath. 

This method may be beneficial for the optimization of operating hybrid processors 
of quantum information comprised of different modules [77, 78, 105]: superconducting 
qubits coupled via a microwave resonator to ensembles of ultra-cold atoms or NV- 
center spins. Hybrid processors may profit from the advantages and make up for 
the shortcomings of the individual modules [105]. Specifically, the superconducting 
qubits are fast but vulnerable to decoherence. The outcome of their operations should 
be controllably transferred to collective “quiet” (decoherence-resilient) states of the 
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atoms that are much better suited for long-term shelving (storage) of this quantum 
information (QI). The overall fidelity of the processor can be improved by dynamical 
control that optimizes this QI transfer from the noisy to the quiet (storage) module. 
Remarkably, for a given energy of the transfer pulse, the shortest pulse is by no means 
optimal [80]! 

A related method can significantly improve the performance of quantum memories 
based on spectrally inhomogeneous spin ensembles [81]. This method preselects an 
optimal portion of the ensemble by appropriate microwave pulse designs. 


3. Bath-induced entanglement in open systems 


Environment effects generally hamper or completely destroy the “quantumness” of 
any complex device. Particularly fragile against environment effects is quantum 
entanglement (QE) in multipartite systems. This fragility may disable quantum 
information processing and other forthcoming quantum technologies [2, 3, 4, 5, 6, 
7, 8, 105]: interferometry, metrology and lithography. Commonly, the fragility of QE 
rapidly mounts with the number of entangled particles and the temperature of the 
environment (thermal “bath”). This QE fragility has been the standard resolution 
of the Schroedinger-cat paradox [14, 106]: the environment has been assumed to 
preclude macrosystem entanglement. But is it inevitable that Schroedinger cats die of 
decoherence (as commonly believed [13, 14])? Or, conversely, can a cat be both dead 
and alive in a thermal bath? 

We shed light on these fundamental issues within the simple model of N spin-^ 
non-interacting particles that identically couple to a thermal oscillator-bath via the 
z-component of their Pauli operators. A single spin in such a model undergoes bath- 
induced pure dephasing [13, 14, 107]. Yet, strikingly [108, 109], an initial product 
state of N z-polarized spins can spontaneously evolve via such coupling to the bath, 
into a Schroedinger-cat state, also known as a macroscopic quantum superposition 
(MQS) or GHZ state [106], nearly deterministically (Eig. 4) 


Ifr) +e®5 |J|)\ /(-frl +e 


V2 


72 


+ (1 - rips, P - 1 


( 12 ) 


GHZ 


with \'([) = Itt •• t)^ IJI) = lii •• i) with only a small probability 1 — p of evolving 
into an incoherent state of the N spins. 

This dynamics of the collective spin along z, is driven by the 

Hamiltonian 

H = looLz + + Lz Vk{bk + 5^) 

k (13) 

"-V-" 

Hi=LzB: collective coupling to bath 


where ooqLz stands for the collective energy (without the bath), bk and respectively 
annihilate and create bath quanta in modes labeled by k, with frequencies uj^ and their 
coupling constants to Lz are denoted by rj^. The evolution of the combined system- 
bath state is exactly soluble by means of the unitary evolution operator [108, 109] 


—i uJotLz + 


fml 


+ Lzt 'y^((^kbk + 


k 


U (t) = exp 


collective Lamb shift- 


(14) 



Thermal Baths as Quantum Resources: More Friends than Foes? 


9 


Noninteracting 
spin/atomic ensemble 



Figure 4. Schematic view of a product-state spin-polarized ensemble (left) that 
spontaneously evolves in the bath into an entangled MQS or GHZ (Schroedinger- 
cat) state at a particular time, as a result of bath-induced entanglement. 



Here the bath-induced nonlinear term is a collective Lamb shift (whose time- 
dependence is described by the bath-dependent function /(t)). This term does not 
affect the bath and only entangles the spins. The last term (wherein Lz is coupled 
to a linear combination of (6|.) with amplitudes a^) gives rise to system-bath 
entanglement which, upon tracing out the bath, decoheres the spin-system state. 

Clearly, the ability to entangle the spins via the bath-induced quadratic Lamb 
shift requires the suppression of the decoherence-inducing term (linear in Lz) in Eq. 
(14). This suppression is achievable by control, consisting of periodic phase flips that 
tend to average out the linear (odd-symmetry) term but leave intact the quadratic 
(even-symmetry) Lamb-shift term. 

Further insight into the competing effects of bath-induced entanglement and 
decoherence can be obtained from a detailed consideration of a realistic model: two- 
atom dispersive coupling to a common cavity bath (Fig. 5a), described by the 
interaction Hamiltonian 


Hi= Yi = Y Y^^ lb) ihlih + bi) (15) 

j=A,B j=A,B k ^ 

Here the energy shift of state |1) in atom j is caused by the combined effect of an off- 
resonant classical field (with Rabi frequency Qj and detuning Aj) and the quantized 
cavity field (with coupling strength g^j) (Fig. 5a). The cross-coupling of atoms 
j = A^B via virtual quanta exchange in the cavity is the source of their collective 
Lamb shift. This cross-coupling is chosen not to depend on the interatomic distance 
under the assumption of identical couplings of both atoms to all cavity modes: 


r]k = 


^3 9k, 3 

^3 


(16) 


which is the case for atoms located at symmetric positions in the cavity. Then the 
foregoing analysis yields the real-quanta exchange rate between the atoms that causes 
decoherence 

^A(B) = 27tGt{(^ = 0) (17) 

where the coupling spectrum of the cavity-bath at temperature T is sampled at cj = 0. 
This decoherence rate competes against the collective Lamb shift 



U3 


( 18 ) 
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Figure 5. (a) Schematic view of bath-induced virtual quanta (wiggly arrows) 
exchanged between atoms A and B, in the presence of off-resonant fields (solid 
arrows). The net result is a collective interatomic energy (Lamb) shift, (b) 
An Ohmic bath spectrum allows for a collective Lamb shift associated with the 
integral over all bath frequencies a;, from 0 to c^c, and thus dominates over the 
decoherence rate associated with the bath at cj ~ 0. 


This two-atom Lamb shift is given by the principal-value part of the integral over 
the entire coupling spectrum, which, remarkably, is taken to be at zero temperature, 
T = 0, regardless of the actual bath temperature. 

The desired dominance of the collective Lamb shift due to virtual quanta exchange 
over decoherence due to real quanta exchange, e.g., in an Ohmic bath, holds if 

Iab > ^A(B) ^ Gt{0). (19) 

Namely, the upper cutoff frequency far exceeds the zero-frequency coupling rate, which 
is typically the case (Fig.5b). 

We thus arrive at the following paradigms: (i) QE in large multipartite systems 
may naturally (spontaneously) arise (albeit over limited time) when the system is 
embedded in commonly encountered thermal environments (baths). This QE may 
yield the spontaneous formation of Schroedinger-cat states (MQS). (ii) QE control 
may actually take advantage of the coupling to the environment rather than try to 
eliminate it, z.e., it should enhance the “helpful” coupling, leading to virtual quanta 
exchange, and suppress the “harmful” exchange of real quanta via the bath. 

Such natural, yet unitary, evolution within thermal baths of the system to 
a highly-nonclassical MQS state is a universal effect which we dub bath-induced 
entanglement (BIE). Whereas, as a rule, the interaction of quantum system with 
a thermal bath gives rise to decoherence, BIE arises from nonresonant (virtual) 
interactions between particles via the bath: nonlinear frequency pulling. This is a 
generalization of effects that have previously been studied for multi-ion coupling to 
single-mode phonons [110]. 

A complementary (orthogonal) approach taken by other groups is to realize 
certain entangled states by engineering the incoherent (nonunitary) dissipation of 
quantum systems into a Markovian (spectrally flat) bath [111, 112]. By contrast, 
the coherent, bath-induced evolution discussed above crucially depends on having a 
non- flat bath spectrum. 

Large ensembles of two-level atoms as considered above may be isomorphic to 
spin systems with large-spin eigenstates. The interaction of such ensembles with a 
common light field may lead to their entanglement [113]. Here, instead, we rely on the 
spectra of commonly encountered baths to drive the spin ensemble into an entangled 
state, via effectively nonlinear dynamics. 

At the same time, we must be concerned with the protection of such bath- 
induced entangled states from the disentangling effects of other baths that constitute 
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Figure 6. Phot on-induced interaction between identical two-level atoms. 
Atom 1, initially in its upper state |ei), emits a quantum at mode k and 
corresponding frequency ujk while making a transition to its lower state 
1 ^ 1 ). Atom 2, initially in its lower state, becomes excited upon absorbing 
the quantum. When the exchanged photon is real, i.e. for Uk — oJa 
where uja is the frequency of the atomic transition, the interaction gives 
rise to dissipative, and hence probabilistic, cooperative-emission. The 
summation over all other virtual-photon-mediated processes, i.e. over all 
transition amplitudes for which ujk ^ yields the quantum-mechanically 
coherent and hence deterministic exchange process of resonant dipole-dipole 
interaction (RDDI). 


their environment. Such protection presents a challenge: how to optimally control 
multiqubit entangled states? Our ability to face this challenge relies on our universal 
approach to multipartite decoherence control [43, 44, 45, 74, 75, 76] (see above). 

4. Long-range bath-induced dispersive interactions 

As argued above, the key to BIE is virtual quanta exchange via the bath. The BIE 
processes considered in Sec. 3 were restricted to identical coupling of all the atoms 
to the bath modes, and hence their collective Lamb shift is distance-independent. 
However, in general this is not the case: the system-bath couplings in the interaction 
Hamiltonian depend on the positions of the individual atoms via the spatial mode 
functions of the bath modes. In free space the mode functions of the photonic 
bath are 3d plane waves, giving rise to real and virtual quanta exchange which 
both decay with interatomic separations r and correspond to Dicke-like cooperative 
emission/absorption and to cooperative Lamb shifts {i.e. resonant dipole-dipole 
interaction - RDDI), respectively [114, 115]. Whereas for interatomic separations r 
longer than the resonant atomic wavelength the real- and virtual-photon processes are 
comparable (scaling as 1/r), in the near-held zone, i.e. for small r, the RDDI retains 
the familiar dipole-dipole scaling as 1/r^, and can greatly exceed cooperative decay. 
Therefore, only in the near-held zone can free-space RDDI lead to predominantly 
deterministic BIE. On the contrary, RDDI-induced entanglement is never deterministic 
at separations beyond the emission wavelength, where incoherent absorption and 
emission render it probabilistic (Eig. 6). 

4 . 1 . Long-range deterministic entanglement via RDDI 

There is a remedy to this state of ahairs that would render BIE via RDDI nearly 
deterministic in the far zone, i.e. for atoms separated by many wavelengths. This 
remedy is based on bath engineering: shaping photon modes at will by changing 
the geometry of the bath. The idea is to consider optical waveguides, such as a 
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Figure 7. Long-range RDDI in a fiber-grating, (a) Atoms (black dots) 
coupled to a fiber with a modulated refractive index (alternating blue 
and black colors) and a grating period of length A. (b) Illustration of 
the dispersion relation of the fiber-grating cj(/c), k being the photon-mode 
wavenumber on the fiber axis: a gap at frequency u and /c = tt/A is opened 
up, with the upper bandedge at the upper cutoff frequency ujco- (c) The 
density of longitudinal photon modes dkjdijj (inverse of group velocity Vg) 
vanishes in the gap and diverges at the bandedge ujco- (d) RDDI mediated 
by the photon modes from (c): the excited-state population of atom 1 is 
plotted using a non-perturbative theory that goes beyond that of Eq. (21) 
[68]. This illustration is plotted for the D1 transition in Rb87 atoms and 
for an inter-atomic distance oi z ^ 16/im (see [68] for more details). The 
population, initially unity, slightly decreases to 0.9663 and then oscillates 
periodically between 0.9663 and 0, similarly to the prediction of Eq. (21). 
This supports a long-distance entanglement generation with concurrence 
C ~ 0.9663 between the atoms at a distance of roughly 20 atomic resonant 
wavelengths, following an interaction duration of about t ~ 1.8ns [68]. 
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Figure 8. Dispersive forces between atoms, (a) Laser-induced forces: an 
off-resonant laser field El illuminates atoms 1 and 2. The scattered field 
between them, proportional to the RDDI strength A, establishes their 
interaction, (b) Vacuum-induced forces (van der Waals and Casimir): here 
the laser is replaced by the electric field of the vacuum fluctuations Eo{z) 
(see text). 


Bragg grating, where the group velocity of guided photonic modes vanishes at the 
cutoff (band-edge) frequency, giving rise to giant enhancement of the mode density 
[65, 68, 116, 69, 70] (Fig. 7a,b,c). Then, considering atomic resonance frequencies 
within the bandgap but very close to its bandedge (cutoff frequency), two consequences 
emerge: 1) The atoms do not exchange real (resonant) quanta (cooperative decay) 
along the waveguide due to the vanishing photon density of states at the atomic 
resonance, thus eliminating their probabilistic, dissipative, interaction. 2) The atoms 
do however, exchange virtual (nonresonant) quanta via RDDI, mediated by all allowed 
waveguide modes (see caption of Fig. 6). Furthermore, the resulting RDDI exhibits a 
strongly enhanced interaction rate (energy) A and effective range scaling as 


fs 




\/l — ^aj^co 


\/l - Wa/i 


Wr, 


( 20 ) 


The above expressions reveal a RDDI whose strength A is much larger than the 
free-space spontaneous emission rate Fj^, and much longer-range than the atomic 
resonance wavelength A^, as the atomic resonance frequency approaches the cutoff 
(bandedge) frequency c^co- Since spontaneous emission, being inhibited into the guided 
modes, remains at the free-space value whereas RDDI now has a giant value, we may 
describe the interatomic exchange of a photon by the effective Hamiltonian that affects 
two-atom entanglement nearly-deterministically, Le., with high fidelity (Fig. 7d), 


ije// = hA y: V L 


( 21 ) 


where {aj ) are the excitation (deexcitation) Pauli operators of the respective 
atoms. 


4 . 2 . Long-range laser-induced forces 

The same principle that allows for the establishment of coherent BIE via RDDI while 
suppressing the incoherent and dissipative process of emission, can be used to create 
long-range conservative forces between atoms. When atoms are illuminated by an 
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Figure 9. Atoms (black dots) coupled to a transmission-line (TL): a TL is 
typically comprised of two distinct conductors, (a) A coaxial TL comprised 
of two concentric conducting cylinders, (b) A coplanar waveguide may be 
seen as an open version of the coaxial line. It is comprised of a central 
conductor and two other grounded conductors (at a distance a from the 
center). 


off-resonant laser, their motion is affected by two processes in analogy to RDDI and 
spontaneous emission: an inter-atomic conservative force, and a diffusive motion due 
to the scattering, respectively. The direct relation to RDDI and spontaneous emission 
can be seen from the following picture (Fig. 8a): the off-resonant laser virtually excites 
the atoms, which, once excited, can either interact coherently via RDDI, leading to 
a distance-dependent cooperative energy shift and hence a force, or emit photons, 
leading to scattering [117]. 

For atoms coupled to a waveguide with a bandgap spectrum as in Fig. 7, and 
illuminated by an off-resonant laser, the resulting laser-induced force follows the same 
RDDI strength and range as in Eq. (20), while the scattering and hence the diffusion, 
is suppressed. Therefore, the dynamics of the motion of atoms in such a system 
are predominantly affected by an extremely long-range conservative force. Such a 
configuration opens the door for the realization and control of many atoms coupled by 
long-range forces, that are expected to exhibit unique thermodynamic features, such as 
inequivalence of statistical ensembles and anomalously slow relaxation to equilibrium 
[62]. 

4.3. Long-range vacuum-induced forces 

Even more dramatic, giant, enhancement is achievable via the control of the bath- 
geometry, for dipolar forces induced by the electromagnetic vacuum, namely, the 
Casimir and van der Waals (vdW) forces. The idea is to consider atoms coupled 
to an electric transmission line (TL), such as a coaxial cable or coplanar waveguide 
[61] (Eig. 9), which support the propagation of quasi-ld transverse electromagnetic 
(TEM) modes. Then, virtual excitations (photons) of these extended modes can 
mediate much stronger and longer-range Casimir and vdW forces than in free-space. 

The unique feature of the fundamental TEM modes is their dispersion-free and 
diffract ion-free Id propagation, revealed by the /c—dependence of their frequency Uk 
and spatial mode function i4/c(r). 
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(a) UlUf, (b) 



Figure 10. TEM-mediated vacuum energy U between a pair of atoms at a 
distance z, coupled to a TL with transverse dimension a = 10“^Ae, where 
Ae is a characteristic dipole-transition wavelength (typical of coplanar- 
waveguide circuit-QED setups), compared to its free-space counterpart 
Ilfs’- (a) near-zone limit z ^ Ae (vdW regime); (b) far-zone limit z ^ Ae 
(Casimir-Polder regime). 


LOk = \k\c, Ufe(r) 


^ikz 

\/A{x,y)L' 


( 22 ) 


k being the wavenumber in the longitudinal waveguide direction 2 : and A the mode 
area. 

The contribution of the TEM modes to the Casimir and vdW potentials can be 
evaluated by fourth order perturbation theory that yields the energy shift of two atoms 
in their lowest (ground) states coupled to the vacuum of the modes in Eq. (22) [118]. 
An alternative approach consist in recalling the expression for the energy of an electric 
dipole, U = —{l/2)aE^^ where a is the dipolar polarizability and E is the field at the 
position of the dipole [119, 120]. Considering the energy of, e.g., atom 2, the electric 
field at its position Z 2 along the transmission line, is comprised of the ordinary vacuum 
fluctuations, £^ 0 (^ 2 ), and those scattered by atom 1 and subsequently arriving at atom 
2, Esc{z 2 )- To lowest order in the scattering, this scattered field is found by the Id 
propagation equation of the TEM modes, driven by the polarization at the location 
of atom 1, Eq{zi) (Fig. 8b) 


{dl + k^)Esc,k{z) = -iioi^lai{LOk)Eo,k{zi)S{z - Zi)/A, (23) 

where k is the wavenumber of the field fluctuations. The solution of the above equation 
yields Esc,k{z 2 ) in terms of ai{uJk)Eo^k{zi)‘ On the other hand, the interaction energy 
between the atoms 1 and 2, deduced from the dipolar energy U of atom 2, is related 
to the scattered-field and hence, to lowest order in the scattering, is given by, 

[/i 2 cx - E] a 2 (wfc) [Esc,k{^ 2 )^o,k{z 2 ) + h.c.] OC - ^ a 2 {LOk)ai{uk) [Eo,k{zi)Eo,k{z 2 ) + h.c.]. 

k k 

(24) 

Finally, treating the vacuum fluctuations E^o,/c(^) ns a quantum field operator, we 
average Eq. (24) with respect to the vacuum state and obtain the vacuum interaction 
energy between the atoms, mediated by the dominant TEM mode of the transmission 
line. 

This enhanced interaction energy may be analytically evaluated in the form of a 
hypergeometric function E of the interatomic distance z = jzi — Z2l scaled to a typical 








Thermal Baths as Quantum Resources: More Friends than Foes? 


16 



Figure 11. Nonadditivity of dispersive interactions, (a) The total energy of 
3 atoms is given by the sum of their pairwise interactions Q-dded 

to their 3-body interaction U123, which is obtained by considering multiple 
scattering events involving all atoms. When U123 is negligible with respect 
to the pairwise summation part, the total energy is called additive, (b) A 
gas of atoms coupled to a transmission line, where nonadditive effects may 
be observed. 


dipolar-transition wavelength Ae- In the near zone, i.e. for ^ much shorter than this 
wavelength, we obtain a modified vdW interaction as compared to free space [61] 


F( 2 :) TT + 167r^ In , (25) 

that falls off very differently with 2 ) than the inverse 6th power, U oc 1 / 2 )^, that 
characterizes the interaction in free space. By contrast, for far-zone distances (well 
beyond Ae), we find [61] 

as compared to the inverse 7th power falloff in free space, U oc Ij. The resulting 
enhancement can be enormous, as seen in Fig. 10. 

An important outcome of the enhancement of Casimir forces in such a TL 
structure may be the onset of nonadditivity of atom-atom interactions [119, 120, 121] 
(Fig. 11a): already at rather low gas densities we expect that 3-body interactions 
1/123 niay become comparable to the strength of their usual pairwise counterpart I/ 12 , 
because of the extension of the mediating photon modes in a TL over long distances. 
This means that the vacuum energy of a many-atom system may not be represented 
by the sum of its pairwise interactions I/ 12 . More specifically, the ratio of 3-body to 
pairwise interaction energies in free space (3d) scales as [119, 120] 


^123 


a 



(27) 


Therefore, the inverse of the polarizability 1/a sets the scale for a typical density 
^ 1/z^ where this ratio is large and nonadditivity is important. For the Id TEM- 
mediated case, as in a TL, we find in the far-zone regime 


Ui23 


a 


oc 


Id 


(28) 
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Figure 12. Cartoon of a quantum system whose thermal equilibrium with a 
bath (top) is changed towards progressive heating by highly frequent quantum 
non demolition (QND) measurements of its energy in the Zeno regime (bottom 
left) or towards progressive cooling by less frequent QND measurements in the 
anti-Zeno regime (bottom right) [48]. 


where a ^ ^/A is transverse dimension of the TL. Then, for the typical case of z ^ a, 
where the TEM modes are dominant and Id behavior prevails, nonaddivity is expected 
to become important at densities much lower than in free space. Among the possible 
consequences of this nonadditivity are drastic modifications of the effective dielectric 
response and the heat capacity of gases coupled to such structures (Fig. 11b). 

5. Thermodynamic control via quantum Zeno Sz anti-Zeno effects 

It is clearly desirable to cool down or purify a qubit at the fastest rate possible to make 
it suitable for tasks of quantum information processing. The standard, straightforward 
way of cooling a system such as a qubit is by equilibrating this system with a cold 
bath. But can one cool qubits faster than their equilibration time? 

Another issue concerning cooling is that it often involves transitions between 
the qubit levels and other auxiliary levels. But what if such auxiliary levels are not 
available? We have shown that these two obstacles may be overcome by exerting 
highly frequent perturbations on a qubit, such as phase shift, or measurements [48] at 
intervals much shorter than the memory (correlation) time of the bath to which the 
qubit couples, and well within the equilibration times [49, 50]. In such a scenario, a 
non-Markovian treatment of the purification must be adopted [48, 49, 50, 107]. 

Specifically, we have experimentally and theoretically demonstrated (in 
collaboration with L. Frydman’s group) the possibility of purifying a qubit coupled to 
a spin bath, by means of repeated noise-induced dephasing that mimics the effect of a 
non-selective (unread) measurements. We have shown (Fig. 13) that the qubit may be 
cooled down to a predetermined temperature that may be much lower than that of the 
bath by means of a suitable controlled dephasing rate, that conforms to the condition 
for the anti-Zeno effect (AZE). By contrast, a dephasing rate that corresponds to 
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Figure 13. Experimental verification [107] of Zeno heating for QND polarization 
measurements (at intervals r = 0.2 ms) and anti-Zeno cooling (at intervals 
T = 1ms), as in the cartoon above (Fig. 12) . Experimental setup: quantum 
system embodied by the nuclear spin of a carbon atom C is in contact 
with a bath embodied by the nuclear spins of 3 protons. Their off-resonant 
frequency mismatch is changed by the measurements, causing polarization 
decrease (heating) or increase (cooling) of C depending on the interval r. 


a quantum Zeno effect (QZE) leads to heating of the qubit [48, 49, 50, 107]. The 
qubit may exist in the state of predetermined purity to which it was driven by the 
measurements or dephasing as long as the entropy of the bath remains constant. A 
violation of this (Bonn) approximation will render the qubit as well as the bath fully 
mixed, i.e. will thermalize them to infinite temperature [49]. 

6. Heat-machine design by system-bath control: quantum 
thermodynamic bounds 

6.1. Work-information tradeoff in the non-Markovian domain 

Open-system manipulations must be optimized within thermodynamic bounds, 
concerning entropy, work and heat production. However, when these manipulations 
are faster than the bath memory time, so that the Markovian approximation does 
not hold, we must revisit and better understand these bounds. Part of the reason is 
that correlations between the system and the bath, which are ignored in standard 
thermodynamics, may play an important role on non-Markovian time scales. In 
particular, as discussed below, they may invalidate the bound set by Szilard, and more 
quantitatively by Landauer, on the tradeoff between information and work [122, 123]. 

System-bath correlation effects. All existing treatments of heat engines 
are based on the assumption that the working-medium (system) is autonomous: its 
evolution (described by a Lindblad-type master equation for the bath-averaged system- 
state ps{t)) suffices for a thermodynamic analysis of an engine [124, 125] driven by 
Hamiltonian Hs{t). Under this standard assumption, the following expression for the 
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Heat will amplify the oscillations 
and produce work 


The World's Smallest Refrigerator 




Push extremely frequently! 


The oscillator will cool down the cold bath 


Figure 14. Cartoon of the “world’s simplest and smallest” universal heat machine 
acting as heat engine (left) or refrigerator (right). It is comprised of a qubit that 
constantly interacts with two baths at different temperatures. Concurrently, it 
is driven by an oscillator that off resonantly modulates (by periodic Stark shifts) 
the qubit resonance frequency. The modulation rate determines whether the qubit 
will act as a heat engines or refrigerator [71]. 


work, VF, is expected to hold [124] for a closed cycle: 

W = ^ tr{psHs}dt. (29) 

The convention is that the work W is negative if it is performed by the system on 
the external piston. According to Lindblad’s H-theorem, or the Kelvin formulation 
of the second law, IT > 0, i.e. the system cannot do work on the piston in a single¬ 
bath setup. Yet, strikingly, according to our results [126] if the cycle is faster than 
the bath memory time, the system may do net work on the piston! To resolve the 
paradox, we contend that, contrary to the standard assumption [124, 125], it is wrong 
to assume that the system is autonomous in the quantum non-Markovian domain: 
the correlations of the system with the bath are then crucial! Accordingly, we show 
that the second law is upheld if we allow for the energetic and entropic cost incurred 
upon decorrelating the entangled system-bath state that exists in thermal equilibrium 
by the measurement that triggers each cycle. The correct description must account 
for the total work during the cycle, evaluated by considering the total state ptot and 
the Hamiltonian Hfot (of tho system and the bath combined) [126]: 

Wtot = (j) tr{ptotHtot}dt. (30) 


The non-negativity of the work Wtot ^ 0, under a closed-cycle {unitary) evolution 
of the Hamiltonian {Htot), in keeping with the second law, is ensured by the inequality 

Wtot = WfAE>0, (31) 

where AE is the measurement cost of the system-bath decorrelation. Equation (30) 
may still allow the system to do net work during the cycle (IT < 0), but it should 
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Figure 15. Components of the universal heat machine depicted in Fig. 14: 
(1) working fluid (qubit system); (2) two baths at different (hot and cold) 
temperatures, permanently coupled to the qubit (via weak coupling); (3) a piston 
(external held oscillator) that periodically modulates the qubit (level distance), 
i.e. uo{t) about the qubit frequency cuq and with modulation frequency (rate) fl, 
and thereby extracts work (at expense of the hot bath) in the heat-engine regime, 
or provides work (in order to cool down the cold bath) in the refrigerator regime, 
depending on [71]. 


compensate for this work by the energy cost /S.E of the system-bath state preparation. 
This cost comes about from changing the mean system-bath correlation energy from its 
negative equilibrium value [48, 49, 50] to zero (or positive value) after the preparation 
(e.g., via a brief measurement or phase-flip). 

6.2. Minimal quantum heat machines 

One of our main targets is the strive to realize the minimal and simplest thermal 
machines in the quantum domain. These may be conceived as (Fig. 14) quantized 
(harmonic-oscillator) “piston” that couples to a single qubit acting as either a quantum 
heat engine (QHE) or quantum refrigerator (QR) on spectrally non-flat baths. 
Floquet analysis of periodically-driven open systems is used to treat their steady- 
state thermodynamics [71, 127]. This formalism aims to separate those distinctly non- 
Markovian (and non- rotating-wave) effects that may cause anomalous thermodynamic 
phenomena on short-time scales [48, 49, 50, 107] from steady-state thermodynamics. 

The simplest variant to be used as a model for the minimal quantum thermal 
machine is a (frequency-modulated) qubit Hamiltonian 

1 27r 

Hit) = -ujit)a^, Lvit + = i^it), (32) 

where uj{t) is the periodic modulation about the qubit resonance cjq with frequency 
Q. Its coupling to two heat baths is given by the interaction Hamiltonian 

Hir,t=a^0iB^ + B^). (33) 

where the bath operators {B^, B^) are respectively associated with hot (T^) and cold 
(Tc) temperatures, as well as with distinct spectra. Using the period-averaged (coarse¬ 
grained) Floquet expansion, we can find the conditions for the coarse-grained density 
operator ps of the qubit to be a steady state of a Lindblad super-operator, 

£P5 = 0, £ = ^£“, (34) 

q 

where q = ±1,±2,... is the harmonic index, a = c, h is the bath index, and 
Cqps describes the steady-state with the qubit resonance c^o shifted by gU, U 
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Figure 16. Separation of the hot and cold bath spectra required for the universal 
heat machine in figs. 14-15. The hot-bath spectrum G^{uj) (red) is that of black 
body radiation, rising with mode density as whereas the cold bath spectrum 
(blue, Lorentzian or Ohmic) extends up to oucut- The periodic modulation 
(TT-flips) of the qubit creates two sidebands (green): the lower overlaps both bath 
spectra, but the upper only overlaps the hot-bath spectrum [71]. 


being the periodic modulation frequency. This decomposition of the Liouvillian 
effectively replaces each bath by multiple g-harmonic “sub-baths” with differently 
shifted coupling spectra. The merit of this equation is that it employs Lindblad 
(completely-positive) dynamics but still allows for non-flat bath spectra. This steady- 
state expansion can serve to evaluate the heat currents exchanged among the multiple 
harmonic “sub-baths” via the qubit. These currents can be controlled and optimized 
by the modulation and the bath-spectrum engineering. Their signs will determine 
whether the machine functions as QHE or QR, without the need for traditional stroke 
cycles (Carnot, Curzon, Otto). We have studied these QHE and QR models in 
optomechanical [128] and spin-ensemble [129] setups. 

We contend that such a periodically-modulated control qubit coupled to both 
baths is a minimal model for QHE and QR (Eig. 15). It is remarkable that such a 
simple model allows for both QHE and QR actions, by contrast to previous models 
that required 3-level [130, 131] or coupled-qubit [3, 4] configurations. 

Under 7r-flip (periodic) modulation of the qubit there are only two dominant 
Eloquet harmonics at u; = c^o i U. If the cold-bath coupling spectrum G^{uj) has 
an upper cutoff, such that the hot-bath spectra G^{uj = cjo + ^) dominates at high 
energies, but G^(cc; = cjo — U) dominates at low frequencies (Eig. 16), then the heat 
flow from the cold to the hot bath (cold current) Jc is proportional to the product of 
the respective hot- and cold-bath spectra, 

JcOcG'*(wo + fi)G"=(wo-fi). (35) 

Under this condition Jc is positive if 

^/i(^o + ^) < ^c(^o - ^) (36) 

where and Uc are the respective thermal bath occupancies. We then obtain QR 
action, namely, heat removal from the cold bath and its dumping into the hot bath. 
The opposite inequality sign implies work production. Thus we have a single control 
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Figure 17. Performance of the universal heat machine depicted in figs. 14- 
16 (under the assumption that the spectra of the two baths are non-overlaping 
(separated)). The efficiency is plotted as a function of the modulation rate Q [71]. 
The efficiency in the heat engine regime rises with Q up to the Carnot bound at 
Q = Qcrit- At higher Q the machine switches over to the refrigerator regime with 
a coefficient of performance that decreases from the Carnot bound. 


parameter^ the modulation rate, ft: for low rates we have an engine (QHE) and for 
high ones a heat pump (QR) (Fig. 17). However, neither the QHE nor the QR action 
will happen if the bath spectra are inappropriate, so that bath-spectrum engineering 
is crucial. 

6.3. Third-Law Issues 

The limits on cooling power and efficiency of this unconventional QR can serve to 
probe the validity of one of the formulations of the Third Law of Thermodynamics by 
Nernst that prohibits the attainability of zero temperature in finite time. To this end, 
we explore the use of a single driven qubit simultaneously coupled to hot and cold 
spectrally non-flat baths. The possibility of cooling a bath down to arbitrarily low 
temperatures, i.e. cooling rate scaling with temperature, is thus a fundamental issue 
that reflects on the applicability of the Third Law.. As a typical example we assume 
that the QR pumps heat into an infinite hot bath, and out of a cold bath whose heat 
capacity is finite cy < oc, resulting in Tc = Tc{t). Strikingly, we may show that 
arbitrarily low temperature may be reached at finite time (non-exponentially fast) by 
the heat pump, for an appropriate (magnon) cold-bath spectrum thereby challenging 
the Third Law [73]. 

7. Conclusions and Outlook 

Progress in technologies such as quantum information processing (QIP) and quantum 
precision measurements (QPM) or metrology is currently restricted by our ability 
to either minimize the environment effects or actively suppress them by “dynamical 
decoupling”. Based on our theoretical and experimental work, we advocate instead 
taking advantage of the environment (bath) as a resource for quantum technologies, 
provided that we optimize its beneficial effects, preferably by non-unit ary open-system 
manipulations that are less restrictive and more robust than unitary operations. 

To this end, we have identified the following generic tools: 

• Universal dynamical control of open quantum systems, be it coherent or 
projective (nonunitary) that can be optimized (in terms of energy investment) 
for the desired operational task and bath spectrum at hand. Such optimized 
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control may conform to one of two paradigms: the quantum Zeno effect (QZE) or 
the anti-Zeno effect (AZE), depending on the task (Sec. 2). In particular, we have 
introduced the use of dynamical control as a means of maximizing the quantum 
Eisher information (or estimation accuracy) regarding the bath spectrum (Secs. 
2.1-2.2 ). We have also resorted to our universal dynamical control for optimizing 
the tradeoff between the fidelity and duration of quantum information transfer 
via noisy and random media which we model as baths (Sec. 2.3). 

• Bath geometry control by mode confinement to one dimension and spectral- 
cutoff design has been shown (Sec. 4) to allow for drastic increase in the range 
and fidelity of bath-induced atom-atom entanglement (Secs. 3, 4.1) and, even 
more dramatically, the giant enhancement of dispersion (Casimir) forces (Sec. 
4.3). These dispersive mechanisms rely on virtual-quanta exchange via the bath, 
which is enhanced by the engineered bath geometry, as opposed to dissipative 
(real-quanta) exchange, which is suppressed by the chosen bath geometry. 

• Bath-spectra engineering has been shown to be a prerequisite for 
thermodynamic control: AZE employed for high-speed qubit cooling (Sec. 5); 
the intriguing possibility of exceeding the Szilard-Landauer bound by taking 
advantage of system-bath correlations (Sec. 6.1); the operation of simple 
(minimal) quantum heat machines based on a periodically-modulated qubit that 
can attain both high efficiency (near the Carnot bound) and power (Sec. 6.2); 
as well as challenging Nernst’s formulation of the Third Law for the cooling of a 
magnon bath towards the absolute zero (Sec. 6.3). 

We are confident that, however intriguing the above results are, they just barely 
“scratch the surface” insofar as bath-assisted quantum processes are concerned. 
Inevitably, such processes are within the realm of quantum thermodynamics. In order 
to be able to benefit from the quantum control tools discussed above, we should 
revisit the foundations of thermodynamics and reformulate its key concepts and laws 
hy (i) removing the system-bath partition; (ii) exploring coherence and entanglement 
effects on thermodynamic variables and (in) substantiating such fundamental effects 
by studies of different realizations: NV-center spins coupled to a spin bath, spin-boson 
models and boson-boson models in optomechanics. 

The corresponding conceptual goals of such future research are foreseen to be as 
follows: 

a) Use bath engineering, i.e. control of its dimensionality, coupling spectrum and 

quantum state, as a key resource in an effort to push the thermodynamic limits 
of quantum device performance: i) long-time non-Markovian behavior is expected 
for qubits near-resonant with an abrupt spectral cutoff of a quasi ID bath 
previously studied by us [132]. Extensions of these effects to finite-temperature 
baths may be prerequisites to pushing thermodynamics into the hitherto 
unexplored strong-coupling regime where system-bath separability breaks down. 
This regime is expected to allow for (partial) reversibility of the entropy and work 
and thereby alter quantum heat engine (QHE), quantum refrigerator (QR) and 
quantum memory device (QMD) performance, ii) Novel methods of controlling 
system-bath coupling by measurements or phase flips at intervals that violate 
Markovianity, can be developed so as to steer the system-bath dynamics towards 
desired outcomes. 

b) Reexamine the work-effciency Carnot limit derived within the system-bath 

separability (weak-coupling) paradigm: Little is known about the strong-coupling 
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regime in thermodynamics, and we may not rule out that it has surprises in store 
insofar as performance bounds of quantum heat machines are concerned [11], since 
the known definitions of heat currents and power output no longer apply in that 
regime. 

c) Discover quantum-operations speed limits in the thermodynamic limit: the 
rates (speed) of quantum information storage and retrieval, cooling and heat 
engine cycles of quantum systems coupled to thermal baths have unknown 
thermodynamic bounds. To understand these bounds, we should extend our 
previous studies of the third law [73] by discovering the scaling of bath cooling- 
rate with temperature T as T ^0. 

To conclude, thermal baths are promising to be “more friends than foes” for exploiting 
the quantumness of systems that couple to such baths, provided that appropriate 
dynamical control and bath engineering are implemented. 
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